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1.  INTRODUCTION 

Distribution  systems  for  manufactured 
hierarchy  of  warehouses  that  stock  goods  for 
houses  and  to  retail  stores,  at  which  demand 
A multi-facility  distribution  system  such  as 
multi-echelon  inventory  system.  Of  course,  such  systems  are  not  limited 
to  warehouses  and  retail  stores.  For  example,  the  factory  which  produces 
the  goods  could  be  part  of  the  system.  The  factory's  inventories  would 
include  raw  materials  and  work  in  process  at  various  stages  of  completion. 
Another  example  would  be  a repairable  item  system  where  there  are  depot 
and  field  repair  stations.  Other  factors  that  can  make  a distribution 
system  more  complex  are  transshipments,  that  is,  the  redistribution  of 
stock  between  warehouses  on  the  same  level,  repair  facilities  at  some  or 
all  distribution  points  (combination  of  repair  and  consumable  products), 
and  exogenous  demand  at  any  facility  in  the  system,  to  mention  only  a few. 


goods  are  often  composed  of  a 
distribution  to  other  ware- 
for  these  goods  originates, 
this  is  generally  called  a 


Once  the  number  of  installations  and  their  locations  have  been 
fixed,  the  multi-echelon  inventory  problem  is  one  of  finding  the  best 
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inventory  policy  for  each  product  at  each  installation.  There  is  a 

lot  known  about  this  problem  but  there  are  also  many  unanswered  questions. 


The  problem  addressed  in  this  paper  pushes  the  decision  process 
back  one  step  by  assuming  neither  the  number  of  installations  nor  their 
locations  have  been  fixed  for  a multi-product  distribution  system.  Organi- 
zations setting  up  new  distribution  systems  or  re-examining  existing  ones 
have  paid  little  attention  to  this  design  problem.  They  have  generally 


fixed  the  number  of  levels  of  distribution  and  the  number  of  installations 


in  some  arbitrary  manner  and  have  then  attempted  to  determine  optimal 
stockage  levels  for  the  given  system. 


In  Section  2 of  this  paper  we  present  a mathematical  model  for  solving 
this  design  problem;  that  is,  a model  for  finding  the  number  of  echelons, 
number  of  installations,  and  where  products  should  be  stocked  for  a multi- 
product, mu] ti-echelon  inventory  system.  The  model  presented  in  Section  2 
relies  heavily  on  the  ability  to  determine  inventory  policies  for  a given 
multi-echelon  system.  A method  for  finding  these  policies,  using  dynamic 
programming,  is  presented  in  Section  3.  Section  4 generalizes  the  model 
presented  in  Section  2 to  include  capacity  constraints  and  presents  a 0-1 
linear  programming  formulation  of  this  generalized  model.  In  Section  5 we 
present  some,  computational  experience  obtained  in  solving  the  two  design 
models  giver  in  this  paper  and  indicate  areas  of  future  research.  This 
introductory  section  concludes  with  an  illustration  of  the  design  problem 
and  a discussion  of  the  premises  on  which  it  is  based. 


1.1  The  Design  Problem 


Figure  1 shows  a three- installation  arborescent-configuration  system. 
Installation  B (warehouse)  obtains  goods  from  a source  (factory)  and  feeds 


the  two  lower-level  installations,  and  k^  (retail  stores),  where 


exogenous  demands,  and  , occur. 
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• 


Source 


Figure  1.  A three- installation,  arborescent- 
configuration  system. 


Figure  2 gives  two  alternative  designs  for  this  system.  In  Figure  2(a) 
installation  B has  been  removed  from  the  system  and  , A2  are  supplied 

directly  from  the  source  of  production.  Figure  2(b)  shows  another  alterna- 
tive design  r.o  the  basic  three- installation  arborescent  system.  A ^ has 

been  removed  and  its  demand  goes  directly  to  B . Since  the  exogenous 
demand  is  still  located,  physically,  near  where  A^  had  been  located,  the 

dashed  line  indicates  that  individual  shipments  are  made  one  at  a time  to 
the  customers  represented  by  • This  can  be  thought  of  as  mail-order 


business. 


Source 


Source 


(b) 


Figure  2.  Alternative  designs  for  the 
system  of  Figure  1. 


The  phrase  echelon  structures  will  be  used  in  referring  to  the 
alternative  designs  for  a multi-echelon  distribution  system.  The  basic 
echelon  structure  will  be  thought  of  as  the  one  containing  the  maximum 
number  of  levels  and  installations  under  consideration.  For  example,  if 
Figure  1 represents  a basic  echelon  structure  then  the  system  has  a maximum 

3 1 

of  2 = 8 possible  alternative  echelon  structures.  In  general,  if  a 

basic  echelon  structure  is  composed  of  p installations,  then  the  number 

of  alternative  echelon  structures  is  2^  . Of  course,  it  might  be  possible 

to  eliminate  immediately  some  of  the  structures  from  consideration  due  to 
the  availability  of  existing  facilities  or  to  geographical,  political, 
economic,  and  other  constraints. 

^This  number  would  increase  if  we  considered  more  situations  than  an 
installation  simply  being  in  the  system  or  not  in  the  system.  For  example, 
a third  possibility  is  to  have  an  installation  carry  a lifetime  supply  of  a 
product. 
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The  multi-product,  multi-echelon  inventory  system  design  problem 
discussed  in  this  paper  is  based  upon  two  premises.  First,  there  is  no 
reason  to  suspect  that  the  minimum  cost  design  will  be  one  in  which  all 
products  use  the  same  echelon  structure.  That  is,  given  a set  of  echelon 
structures  under  consideration,  if  we  were  to  find  the  best  inventory 
policy  for  each  product  stocked  under  each  echelon  structure,  it  would 
very  likely  be  found  that  different  echelon  structures  were  best  for  dif- 
ferent products.  Second,  the  objective  to  minimise  some  measure  of  total 
system  cost  will  usually  be  accomplished  only  if  some  products  do  not  use 
the  echelon  structure  best  for  them.  This  is  because  different  products 
stored  at  the  same  installation  will  share  the  fixed  cost  of  that 
installation. 

To  illustrate  these  premises  we  have  constructed  an  example  problem 
concerned  with  designing  a four-product  inventory  system  in  which  four 
echelon  structures  are  considered  [see  Pinkus  (197 .b) ] . In  this  example, 
which  uses  realistic  inventory  costs,  we  find  that  a different  echelon 
structure  is  best  for  each  product.  Furthermore,  to  minimize  the  total 
cost  of  the  system,  two  of  the  products  cannot  be  stored  in  ways  that 
minimize  their  individual  inventory  costs. 

The  next  section  presents  a mathematical  model  for  solving  the  design 
problem  described  in  this  section. 

2.  ANALYSIS  OF  THE  DESIGN  PROBLEM 

In  this  section  a mathematical  model  for  determining  the  best  design 
for  a multi-product,  multi-echelon  distribution  system  is  presented.  The 
purpose  of  this  model  is  to  find  the  best  echelon  structure  for  each  product, 
bearing  in  mind  that  the  products  are  not  independent  when  it  comes  to 
sharing  installations.  A superposition  of  the  best  echelon  structure  for 
each  product  will  then  result  in  the  best  system  design. 

2.1  Measure  of  Effectiveness 

We  first  consider  what  is  meant  by  "best"  design.  The  measure  of 
effectiveness  to  be  used  will  be  one  that  minimizes  some  measure  of  total 


system  cost,  for  example,  total  expected  cost  per  year  or  total  discounted 
cost  for  the  lifetime  of  the  system.  The  system  cost  will  include 
all  the  real  ( out  of  pocket  ) operating  expenses  of  the  system, 
capital  expenditures  for  building  the  system,  and  those  intangible  factors 
to  which  a cost  can  be  assigned.  The  penalty  cost  for  loss  of  goodwill 
as  a result  of  back  ordering  (not  meeting  demand  when  it  occurs)  is  an 
example  of  such  an  intangible  cost.  Any  intangible  factor  to  which  a cost 
can  be  assigned  can  be  handled  directly  by  the  model  to  be  presented.  Other 
ways  of  handling  intangible  factors  will  be  discussed  at  the  end  of  this 
section. 

For  a given  product  using  a given  echelon  structure,  we  separate 
these  system  costs  into  two  categories.  The  first  includes  all  the  costs 
associated  with  the  inventory  stockage  policy.  We  call  these  the  inventory 
costs  and  they  include  the  costs  of  procurement,  carrying  inventory,  fill- 
ing orders,  and  stockouts.  The  second  category  includes  all  the  costs  asso- 
ciated with  operating  the  installations  of  an  echelon  structure  regardless 
of  the  number  of  products  using  the  installations.  We  call  these  facility 
costs  and  they  include  the  capital  expenditure  for  building  the  installa- 
tions, along  with  a number  of  fixed  costs  associated  with  operating  the 
installations,  for  example,  administrative  expenses,  the  expense  of  renting 
facilities  (if  they  are  not  built),  and  certain  other  fixed  operating  expenses 
for  a given  product  that  do  not  depend  on  the  inventory  policies  used  at 
the  various  installations.  The  facility  costs  are  all  fixed  costs  that  are 
normally  considered  sunk  costs  in  solving  a multi-echelon  inventory  problem, 
because  the  echelon  structure  has  been  set. 

2.2  Mathematical  Model  of 
the  Design  Problem 

Let  a = the  inventory  cost  of  product  j 
using  echelon  structure  i , 

b^  = the  facility  cost  of  installation  k , 

= the  facility  cost  of  installation  k 
when  echelon  structure  i is  used. 
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where 


h,,  * 
lk 


b^  , if  installation  k 
echelon  structure  i , 

0 , otherwise. 


is  included  in 


The  costs  a.  ^ and  are  arrayed  in  the  matrix  [a.  . :b.v]  of  Figure  3 


ij*  ik 


Products  (j) 


Installations  (k) 


Echelon 


Structures 


(i) 


1 

2 

• • • 

n 

1 

2 

• • • 

P 

11 

a12 

• • • 

3ln 

' bll 

i 

b12 

• • • 

blp 

21 

a22 

• • • 

a2n 

H 

CM 

b22 

• • • 

b2p 

• 

1 

1 • 

• 

• • 

• 

• • 

1 

ml 

am2 

• • • 

a 

mn 

1 b , 
i mj. 

bm2 

• • • 

b 

mp 

Figure  3.  Array  of  inventory  costs  and  facility  costs. 


If  the  facility  costs  b^  were  all  zero,  all  products  could  be  treated 

independently  and  the  solution  to  this  problem  would  be  obvious.  We  would 
select  the  minimum  element  in  each  column  of  matrix  [a„]  » that  is,  the 

minimum  inventory  cost  for  each  product.  The  echelon  structure  i 
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associated  with  this  minimum  cost  for  product  j would  he  the  echelon 
structure  used  by  product  j . Unfortunately,  the  b_.,  are  not  all  zero. 

IK. 

The  design  problem  is  similar  to  the  assignment  problem  of  linear 
programming  in  that  it  is  desired  to  assign  products  to  echelon  structures 
so  that  minimum  cost  over  the  is  achieved.  Unlike  the  assignment 

problem,  however,  multi-assignments  in  any  row  are  permitted  and  the  assign- 
ment of  products  to  echelon  structures  can  create  additional  costs  (facility 
cost),  depending  upon  which  installations  are  required  in  the  echelon  struc- 
tures chosen.  A further  complication  occurs  because  each  facility  cost, 
b^  > is  a fixed  cost,  regardless  of  the  number  of  products  using  installa- 
tion k . 

We  now  define  decision  variables  and  give  a mathematical  formulation 
for  the  design  problem. 

Let  x„  = 1 if  product  j uses  echelon  structure  i , 

=0  if  product  j does  not  use  echelon  structure  i . 

Then  the  problem  is  to  find  the  matrix  X = [x..l  which  minimizes 

ij 


m n p 

F(X,'  = Z Z a,,x_,,  + Z u 


i=l  j=l 


ij  ij 


k=l 


2 (bik  Z xij) 
i=l  j=l  J 


(1) 


subject  to: 


m 

Z 

i=l 


xij  = 1,  j=l,2, . 


x.  . = 0 or  1 , for  all  i and  j , 
ij  J 


(2) 


(3) 


where 


u(y)  =0  if  y _<  0 , 
= 1 if  y > 0 . 
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Thus  the  problem  of  designing  a multi-product,  multi-echelon 
distribution  system  has  been  formulated  as  a 0-1  nonlinear  programming 
problem.  The  objective  is  to  find  the  matrix  X*  which  minimizes  F(*)  . 

Real-life  design  problems  often  must  take  into  account  intangible 
factors  to  which  costs  cannot  be  applied.  For  example,  suppose  it  is 
necessary  to  have  a particular  installation  in  the  system  because  of  cer- 
tain trade  agreements,  regardless  of  the  cost  this  might  entail.  This 
situation  may  be  handled  by  considering  only  those  echelon  structures 
which  include  the  particular  installation. 

We  have  solved  the  above  mathematical  model  using  a branch-and-bound 
algorithm.  This  algorithm  is  described  elsewhere  [see  Pinkus,  Gross,  Soland 
(1973)].  However,  one  of  the  weaknesses  of  the  formulation  just  presented 
is  that  it  assumes  that  there  is  no  limit  to  the  storage  space  available  at 
a given  installation.  It  is  one  purpose  of  this  paper  to  overcome  this 
weakness  by  incorporating  space  constraints  in  this  model.  Before  doing 
that,  however,  we  describe  one  way  of  finding  the  inventory  costs,  a„  . 

To  obtain  these  costs  requires  the  solution  of  another  problem — the  multi- 
echelon inventory  problem. 

3.  MULTI- ECHELON  INVENTORY  MODEL 

This  section  is  concerned  with  obtaining  the  value  of  a„  , the 

2 

inventory  cost  to  be  used  in  the  design  model.  This  value,  and  the 
associated  inventory  stockage  policies,  are  arrived  at  by  solving  a multi- 
echelon inventory  problem.  Thus,  for  product  j stocked  under  echelon 
structure  i , it  is  desired  to  find  the  optimal  inventory  policies,  at 
each  installation  of  the  structure,  which  yield  a„  . 

Only  dynamic  solutions  to  the  multi-echelon  inventory  problem  were 
considered.  This  means  that  we  want  to  make  a number  of  inventory  stockage 
decisions  over  time,  not  consider  just  a one-period  or  static  solution. 


To  avoid  confusion  in  this  section,  the  notation  a.,  will  be  used 


ij 


rather  than  the  phrase  inventory  cost,  where  the  particular  echelon  structure 
i or  the  particular  product  j is  immaterial.  The  methods  for  obtaining 
a,,  will  be  applicable  for  all  i and  for  all  j . 


ij 
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At  first  glance,  this  seems  to  be  an  obvious  assumption  for  the  inventory 
problem  so  that  it  will  tie  into  the  problem  of  designing  a system  to  be 
used  for  a long  period  of  time.  However,  it  is  conceivably  optimal  for 
a distribution  system  to  stock  a lifetime  supply  of  a product,  rather 
than  reorder  from  the  source  of  production.  If  one  wanted  to  consider 
such  a situation,  it  would  best  be  handled  by  modifying  the  design  model. 

3.1  Clark's  Approach 

Clark  (1958)  presented  a dynamic  programming  solution  to  the 
problem  of  finding  inventory  policies  for  a multi- installation,  single- 
product inventory  model  with  stochastic  demand.  His  method  for  solving 
multi-echelon  inventory  problems  is  built  on  the  framework  of  the  classical 
approach  to  uncertain  demand,  periodic  review  dynamic  inventory  problems, 
first  presented  by  Arrow,  Harris,  and  Marschak  (1951).  Under  this  approach, 
the  cost  of  the  system  is  represented  as  a function  of  the  inventory  level. 
This  cost  function  includes  ordering,  holding,  and  shortage  costs.  The 
objective  is  to  minimize  the  total  expected  discounted  cost. 

Clark's  solution  was  shown  to  be  optimal  by  Clark  and  Scarf  (1960) 
when  the  installations  were  in  series.  However,  they  pointed  out  that 
Clark's  procedure  was  not  necessarily  optimal  for  the  more  realistic 
arborescent  configuration  of  installations  (see  Figure  1).  Nevertheless, 
they  felt  that  under  many  real  situations  it  gave  an  excellent  approximation 
to  the  solution  of  the  arborescent  case. 

One  of  the  assumptions  made  by  Clark  and  Scarf  was  that  the  cost 
to  an  installation  of  ordering  an  item  from  a higher  installation  in  the 
system  was  linear,  without  any  fixed  cost  of  ordering.  The  only  exception 
to  this  assumption  was  at  the  highest  installation,  where  a fixed  cost  of 
ordering  was  allowed.  In  a subsequent  paper  [Clark  and  Scarf  (1962)]  they 
relaxed  this  assumption  and  were  able  to  give  upper  and  lower  bounds  on 
the  cost  function  for  the  optimal  solution.  This  led  to  an  approximate 
optimal  inventory  policy  for  each  installation,  ased  upon  the  upper  bound 
on  the  optimal  solution. 
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A modification  of  the  Clark  approach  will  be  used  to  determine  the 
values  of  a^  for  the  design  model  of  this  paper.  It  should  be  noted 

that  any  method  for  determining  multi-echelon  inventory  policies  could  be 

used  to  find  the  values  of  a..  . The  selection  of  a method  is,  of  course, 

13 

dependent  upon  the  characteristics  of  the  inventory  situation  for  which  a 
multi-echelon  system  is  being  designed.  The  Clark  approach  has  been 
selected  because  (1)  the  determination  of  the  a.,  is  an  immediate  by- 

13 

product  of  the  calculation  of  the  optimal  inventory  policies,  (2)  the 
dynamic  programming  calculations  are  straightforward,  and  (3)  the  model 
is  rich  enough  to  include  many  multi-echelon  inventory  situations. 

3. 2 Review  of  Classical  Inventory  Models 

Before  describing  Clark's  method  for  finding  optimal  inventory 
policies  for  a single  product  stored  in  a multi-echelon  system,  we  review 
the  classical  dynamic,  periodic  review,  stochastic  demand  inventory  model 
for  a single  item  at  a single  installation.  This  model  forms  the  basis 
for  the  Clark  multi-echelon  model. 

3.2.1  Classical  Inventory  Problem  - A purchasing  decision  is  to  be 
made  at  the  beginning  of  each  of  a number  of  regularly  spaced  periods  of 
time,  for  example,  at  the  beginning  of  each  week.  This  decision  will  be 
based  on  the  level  of  inventory  at  the  time,  ordering,  holding,  and  shortage 
costs  during  the  period,  and  the  effect  the  decision  will  have  on  future 
periods.  Let 

z = purchase  quantity,  z >_  0 

c(z)  = cost  of  purchasing  z units 

X = number  of  time  periods  lag  between  an  order  and 
its  delivery,  the  possible  values  for  X being 
0,1,2, ...  . 
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4>(t)  = probability  density  function  of  demand  during 
a period,  where  demand  is  a continuous  random 
variable  and  is  independent  from  period  to 
period 

u = inventory  on  hand  at  the  end  of  a period,  where 
-°°  < u < 00  (a  negative  value  indicates  that 
demand  occurred  during  the  period  that  could  not 
be  filled) 

h(u)  = holding  cost  charged  on  inventory  on  hand  at  the 
end  of  a period 

p(u)  = shortage  cost  charged  for  failure  to  meet  demand 
during  a period 

x = inventory  on  hand  at  the  beginning  of  the  nth 
period,  before  an  order  is  received,  that  is, 
the  inventory  on  hand  at  the  end  of  the  previous 
period 


y = inventory  on  hand  at  the  beginning  of  the  nth 
period,  immediately  after  an  order  is  received 

L(yn)  = expected  holding  and  shortage  cost  during  the  nth 
period  (hereafter  referred  to  as  the  period  cost). 

/Qn  h(yn-t)<|>(t)dt  + Jy  p(t-yn)<Ht)dt  , yn  _>  0 


r°° 

Jq  p(t-yn)<j>(t)dt  , yn  < ° 


n 1 , 2 , . . . 


If  a delivery  is  to  be  received  during  a given  period  as  a result  of 
an  order  placed  A periods  before,  then  it  is  assumed  this  order  arrives 
at  the  beginning  of  the  period  and  before  the  purchase  decision  for  the 
period  is  made.  Furthermore,  it  is  assumed  the  supplier  carries  an  infinite 
supply  of  the  item,  that  is,  the  supplier  never  back  orders  the  installation. 
The  dynamic  programming  formulation  of  this  problem  is  now  given.  In  this 
formulation  it  is  assumed  that  A = 0 , that  is,  delivery  is  instantaneous, 
and  that  excess  demand  is  back  ordered.  Let 


C (x  ,y  ) = total  expected  discounted  inventory  cost 
for  a problem  lasting  n periods 


1,2, 


a = discount  factor. 
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The  periods  are  numbered  backwards  in  time,  thus,  period  number  one  is 
the  last  period  of  the  problem.  It  is  assumed  that  units  on  hand  at  the 
end  of  the  last  period  have  no  salvage  value. 

Suppose  we  are  at  the  beginning  of  the  nth  period  of  the  problem, 

that  is,  there  are  n periods  of  business  remaining  for  the  installation, 

and  x is  the  inventory  on  hand  before  an  ordering  decision  is  made.  The 
n 

optimal  policy  for  the  nth  period  is  the  policy  which  minimizes  C^Cx^.y^)  * 
The  well-known  dynamic  programming  recursive  relation  for  this  problem  is: 


C (x  ) = min  (c(y-x)  + L(y  ) 
n n . n n n 

y > x 
n—  n 


n = 1,2,...  (4) 


+ a /q  Cn-l<Vt)<Ht)dt} 


where  C^Cx^)  equals  minimum  total  expected  discounted  cost  for  a problem 

lasting  n periods,  n = 1,2,...  . In  this  equation  C (x  ) has  been  broken 

n n 

down  into  three  components:  the  purchasing  cost  for  the  nth  period;  the 
period  cost  for  the  nth  period;  and  the  total  expected  discounted  cost  for 
n-1  periods  of  operation,  assuming  an  optimal  inventory  policy  is  followed 
during  the  last  n-1  periods.  This  recursive  relation  is  used  to  find  the 
optimal  value  of  y^  , which  we  call  . Clearly,  the  desired  inventory 

level  at  the  beginning  of  the  nth  period,  , has  an  effect  on  all  future 
levels  Sj.  , i = 1,2,..., n-1  . 

Two  well-known  results  for  this  model  are  now  given. 


3.2.2  No  Fixed  Cost  of  Ordering  - Assume  the  purchasing  cost  is 
linear  with  no  fixed  cost  of  ordering.  Then 

c (z ) = c • z , z •>  0 . 

Assume  L(yn)  is  convex.  Then  it  has  been  shown  by  Karlin  (1958)  that 

the  optimal  policy  for  an  n period  problem  can  be  characterized  by  a 
sequence  of  critical  numbers  ’ T^e  f°r  c^e  kth  period 

is: 
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if  <_  , order  , 


if  x^  > S^  , order  nothing. 


This  result  can  be  shown  as  follows: 

Vvyk>  " c(Vxk>  + L<yk)  + ° Q ck.1(yk-t)‘Kt)dt  . (5) 

C (y  -t)  is  a convex  function  because  it  is  the  sum  of  convex  functions. 

k“  1 k 


Let 


/CO 

o ck_i(Vt)0(t)dt 


W,  (y,  ) is  the  sum  of  convex  functions  and,  thus,  is  a convex  function, 
k k 


c(yk-xk>  = c * (yk-xk)  • 


Substituting  into  (5) 


Ck(xk>yk)  = C * (yk-Xk}  + Wk(yk}  * 

We  desire  to  find  , the  value  of  which  minimizes  function  . 

Taking  the  partial  derivative  and  setting  equal  to  zero,  we  get 


W.  (S,  ) = -c  . 
k k 


3.2.3  Positive  Fixed  Cost  of  Ordering  - Assume  the  purchasing 
cost  is  linear  with  a fixed  cost  of  ordering  equal  to  K . Then 

c(z)  =c*z+K,z>0, 

= 0 , z = 0 . 

Assume  -*-s  convex.  Scarf  (1960)  has  shown  that  the  optimal  policy 

for  the  kth  period  is  defined  by  a pair  of  critical  numbers,  (Sk»sk)  • 

The  policy  for  the  kth  period  is: 

if  x,  < s,  , order  S - x,  , 
k — k k k 

if  xk  > sk  » order  nothing. 
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3.2.4  Two  Unnecessary  Assumptions  - In  presenting  the  dynamic  program- 
ming formulation  for  the  classical  single-installation  model,  two  unnecessary 
assumptions  have  been  made  for  the  sake  of 

(a)  simplifying  the  following  description  of 
Clark's  approach,  and 

(b)  simplifying  the  computer  program  which  has 
been  written  to  obtain  the  inventory  costs, 


First,  it  was  assumed  that  A = 0 , that  is,  delivery  of  an  order 
is  immediate.  This  is  not  a necessary  assumption.  For  the  case  where 
A > 0,  Karlin  and  Scarf  (1958)  show  that  the  optimal  purchase  quantity  is 
a function  of  the  total  stock  on  hand  plus  on  order,  regardless  of  the 
dates  of  delivery,  assuming  excess  demand  is  backlogged.  They  show  how  to 
reduce  a problem  with  an  order  time  lag  (A>0)  to  one  in  which  no  lag  exists, 
and  thus  the  above  results,  where  A = 0 , are  applicable  in  general.  In 
developing  their  multi-echelon  results,  Clark  and  Scarf  allow  for 
A _>  0 (A=0,1, 2, . . . ) . To  simplify  the  discussion  of  these  results  we 
continue  to  assume  A = 0 . 

Second,  it  has  been  implied  by  the  lack  of  a subscript  that  the 
functions  <Kt),  c(z),  h(u),  and  p(u)  are  stationary.  The  dynamic  pro- 
gramming formulation  for  solving  the  single-installation  problem  can  easily 
handle  a situation  where  these  functions  are  non-stat ionary.  Nevertheless, 
the  assumption  that  they  are  stationary  is  continued.  This  greatly  simpli- 
fies the  amount  of  input  data  and  bookkeeping  necessary  for  the  operation 
of  the  computer  program  that  has  been  written  to  determine  the  inventory 
costs.  Finally,  it  is  noted  that  although  4> ( t ) was  stated  to  be  a con- 
tinuous density  function,  all  the  results  of  this  section  hold  if  (t) 

. 

is  a discrete  probability  function.  The  following  discussion  assumes  the 
continuous  case,  while  the  examples  and  computer  program  utilize  the 
discrete  case. 
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3.3  Multi-echelon,  Series-configuration  Model 

Having  reviewed  two  important  cases  of  the  classical  single-echelon 
model,  we  are  now  in  a position  to  describe  the  Clark  approach  to  a multi- 
echelon stochastic  demand  model.  For  the  present,  only  a multi-echelon 
system  in  which  the  installations  are  arranged  in  series  is  considered. 

Suppose  there  are  N installations,  where  installation  N supplies 
stock  to  installation  N-l  , N-l  supplies  stock  to  N-2,...,  installation 
2 supplies  stock  to  installation  1 . This  series  configuration  is 
pictured  in  Figure  4. 


Source 

[N] 


[N-l] 


vl' 

[2] 


[1] 

tt 

Demand 


Figure  4.  A multi-echelon  system  with  N installations 
in  a series  configuration. 
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The  highest  installation  in  the  series,  N , receives  its  stock  from  the 

3 

source  of  production.  The  notation  and  assumptions  given  at  the  beginning 
of  this  section  still  generally  apply  and  are  not  repeated  here,  with  the 
following  exceptions: 

Let  x1  = inventory  on  hand  at  echelon  i at  the  beginning 

of  the  nth  period  before  an  order  is  received 

y*  = inventory  on  hand  at  echelon  i at  the  beginning 
of  the  nth  period,  immediately  after  an  order  is 
received 

L.Cy1)  = period  cost  of  echelon  i 
1 r. 

ci(z)  = ordering  cost  function  for  echelon  i , z ^ 0 . 

It  is  important  to  note  the  following  distinction  between  an  installation 
and  an  echelon.  The  stock  at  installation  i refers  only  to  the  stock 
physically  at  that  location.  But  when  we  refer  to  the  stock  at  echelon  i , 
we  mean  the  sum  of  all  the  stocks  at  installations  i,  i-l,...,2,l  plus 
all  the  stocks  in  transit  between  installations  i,  i-l,...,2,l  . 

The  formulation  of  the  Clark  model  follows.  It  is  assumed  that: 

1.  Demand,  exogeneous  to  the  system,  occurs  at 
installation  1 only. 

2.  The  purchasing  cost  between  installations  is 
linear  without  a fixed  cost  of  ordering.  This 
cost  can  be  thought  of  as  the  cost  to  transport 
a unit  from  one  installation  to  the  next  instal- 
lation in  the  system.  The  only  exception  to 
this  assumption  is  at  the  highest  installation, 
where  a fixed  cost  of  ordering  is  allowed. 

3.  Demand  in  excess  of  supply  at  any  installation 
is  backlogged. 


3 

The  source  of  production  could  be  included  in  the  series  as  the 
highest  installation.  Then,  its  source  of  stock  would  be  the  supplier  of 
raw  material. 


i 
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4.  The  functions  L^(y^)  are  convex,  i = 1,2,...,N  . 

4 

5.  Delivery  at  each  installation  is  instantaneous. 

6.  The  functions  c(z),  h(u),  p(u)  and  cj)(t)  are  stationary. 

The  model  is  described  for  the  special  case  N = 2 . Let 
~ 1 2 

C (x  ,x  ) = minimum  total  expected  discounted  cost  of 
the  system  at  the  beginning  of  period  n 
before  an  order  is  placed,  n = 1,2,...  . 

Following  the  approach  for  the  single- installation  model,  we  get  the 
recursive  relation: 


: , 1 2.  . , , 1 1,  . ,2  2. 

C (x  ,x  ) = min  (c. (y  -x  ) + c(y-x) 

n n n ^ ^ 2 inn  z n n 

x < y < y 
n — 'n  — n 

2 . 2 
y > x 

n — n 


n — 1 , 2 , • • • 

+ Li (y^  + L?(y^  + a In  Cr>  i(yl-t»  y^-t)4>(t)dt} 

In  l n u n-l  n n 


The  problem  with  this  approach  is  that  in  the  general  case  is  a 

function  of  N variables.  Therefore,  the  recursive  calculations  of 
dynamic  programming  would  be  prohibitively  long,  even  for  a high-speed 

/v 

computer,  if  the  function  Cn  is  left  in  this  form.  Clark  and  Scarf 
have  overcome  this  difficulty.  They  prove  that  the  function 


~ 1 2 N 

C (x  ,x  x ) can  be  decomposed  into  N functions,  each  of  a single 

n n n n 

variable,  that  is, 


; , 1 2 n.  :i,  i.  . :2,  2.  . :n,  n. 

C (x  ,x  , . . . x ) = D (x  ) + D (x  ) + ...  + D (x  ) 
nnn  n nn  nn  nn 


This  permits  the  computation  of  the  optimal  inventory  level  at  each  echelon 
separately  using  the  same  method  that  is  used  for  the  classical  single- 
installation model.  What  is  involved  is  described  by  returning  to  the 
special  case,  N = 2 . 


It  was  pointed  out  earlier  in  this  section  that  this  is  not  a 
necessary  assumption. 
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~ \ 2 ~i  1 2 

C (x  ,x  ) = D (x  ) + D (x  ) , n = 1,2,...  . 
n n n n n n n ’ ’ 


Begin  by  considering  the  lower  echelon  and  assume  there  is  no  limitation 
on  the  stock,  available  to  it  from  the  higher  echelon. 


dV)  = min  {c  (yW")  + L (y*) 
n n . inn  In 


i _ i 
y > x 
n — n 


+ a Iq  D^_1(y^-t)4>(t)dt}  , n = 1,2,. 


(6) 


It  has  been  assumed  that  L^(yn)  is  convex  and  the  purchasing  (transportation) 

cost  is  linear  without  a fixed  cost  of  ordering,  that  is, 

c. (y^-x^)  = c,  • (y^-x^)  . Therefore,  recursive  relation  (6)  is  identical 

J.  n n i n n 

to  Equation  (4),  and  the  optimal  policy  for  this  echelon  is  described  by 

the  sequence  of  critical  numbers  , where  is  the  optimal 

level  of  inventory  for  echelon  1 in  the  ith  period.  This  conclusion  has 
assumed  installation  2 can  satisfy  the  demand  of  installation  1.  This  will 
not  always  be  the  case.  However,  Clark  and  Scarf  show  that  installation  2 
should  satisfy  as  much  of  the  demand  from  installation  1 £ s is  possible 
and  backlog  the  rest.^ 

The  fact  that  installation  2 might  not  be  able  to  satisfy  the  demand 
of  installation  1 suggests  that  in  solving  for  the  optimal  policy  of  eche- 
lon 2,  a shortage  penalty,  in  addition  to  the  penalty  included  in  L ^ , be 

incurred  by  installation  2 when  it  must  backlog  the  demand  that  installation  1 
places  on  it.  This  penalty  is  simply  the  additional  expected  cost  suffered 
at  installation  1 because  installation  2 could  not  satisfy  its  demand.  This 
penalty  is  determined  as  follows. 


It  is  assumed  that  installation  2 will  deliver  at  most  one  shipment 
each  period  to  installation  1 and  that  this  shipment  goes  out  in  the  "split 
second"  after  installation  2 has  made  its  reorder  decision  but  before  its 
stock  is  replenished.  This  sequence  of  decisions  is  necessary  in  order  to 
force  this  example,  with  X = 0 , to  behave  like  a real  problem  with  X > 0 . 
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Consider  the  kth  period  at  echelon  1.  Suppose  < sj  . Then 

k k 

echelon  1 will  order  - xd  stock  from  echelon  2,  and  the  total  cost 

K K. 

for  k periods  at  echelon  1 will  be 

“k<*k>  ■ ci  • <sK> + h(i> + “ r0  . 

2 11 

If  the  stock  level  at  echelon  2 is  such  that  x < (S  -x,  ) , that  is, 

k.  k k 

installation  2 cannot  satisfy  all  the  demand  from  installation  1,  then 

2 1 
echelon  1 will  have  x^  stock  on  hand  after  ordering,  rather  than 

stock.  This  means  that  the  order  received  by  installation  1 is 
2 1 

x^  - x^  . In  other  words,  all  the  available  stock  at  installation  2 was 

shipped  to  installation  1.  Under  this  situation,  the  total  cost  for  k 
periods  at  echelon  1 will  be 


"k(VXk}  = C1  * (VXk)  + Ll(x^  + a /o  Dk-l(xk_t)<f)(t)dt 


Let  A (•)  = the  additional  shortage  penalty  at  echelon  i 
for  not  being  able  to  meet  demand  at  echelon 
i-1  during  the  nth  period,  i = 2,3,...,N  , 
n 1,2,...  . 


The  penalty  A^( • ) is 

Bk<xk’xk>  - Dk(xk> 


<v\>  - ci  • <sH> + Li(xk>  - Li(si> 

+ <*  Iq  [Dk-i(xk_t)  " Dd_1(Sd-t)]4>(t)dt 


= cx  • (x^-Sd)  + L(x^)  - L(Sd) 


+ a/o  - Di_i(S^t)]$(t)dt 


1 


2 2 

W • 
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A more  complete  statement  of  the  function 


is 


2 2 2 1 2 1 
A fx  ) = c • (x^-sb  + L (x  ; - L/S1) 
nn  1 nn  In  lri 


+ aQ  '“i-i'V*)  - . 


for 

x2  < 

sJ 

0 , for 

n — 

X2  > 

n 

s1 

n 

n 

n 1 1 2 } • • • 

n 1 j 2 1 1 * t • 


The  fact  that  A is  a function  of  x alone  is  very  significant.  It 
n n 

means  that  the  recursive  relation  for  the  total  cost  of  n periods  of 

2 

operation  at  echelon  2 is  a function  of  the  single  variable  x . Thus, 

n 

a dynamic  programming  solution  is  feasible. 

The  iptimal  policy  at  echelon  2 is  now  described. 


~2  2 

D (x  ) = min 
n n 2 2 

y > x 

n — n 


{c2<V*n>  + L2(*n>  + An(*n> 


+ <*  !q  D^_1(y^-t)<})(t)dt}  , n = 1,2,. 


(7) 


It  has  been  assumed  that  l^Cy  ) is  convex  anc^  t^ie  purchasing  cost  is 

2 2 2 2 

linear  with  a fixed  cost  of  ordering,  that  is,  c„(y  -x  ) = c„  • (y  -x  ) + K 

2 n n 2 n n 

2 

Since  An  is  the  sum  of  convex  functions,  it  is  a convex  function.  There- 
fore, the  solution  of  (7)  follows  the  solution  given  for  the  classical  model 
with  a fixed  cost  of  ordering,  and  the  optimal  policy  for  the  nth  period 

2 2 

at  echelon  2 is  described  by  the  pair  of  critical  numbers  (S  ,s  ) , where 

n n 

the  superscript  indicates  the  echelon. 

The  procedure  described  here  for  two  echelons  can  be  used  for  any 
number  of  echelons  in  series.  Clark  and  Scarf  show  that  the  policies 
obtained  in  this  way  are  optimal  policies,  and  their  proof  allows  for 

A > 0 . 
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A summary  of  the  important  ideas  behind  the  Clark  approach  follows. 

1.  The  decision  variable  is  the  stock  level  of  an  echelon, 
not  the  stock  level  of  an  installation. 

2.  An  N installation  problem  in  N decision  variables 
becomes  N separate  problems,  each  in  one  decision 
variable.  Each  separate  problem  is  a special  case  of 
the  classical  single- installation  stochastic  demand 
inventory  problem  and  can  be  solved  by  dynamic  programming. 

3.  The  cost  function  for  each  echelon  i (except  the  first) 

includes  the  additional  shortage  penalty  A^  for  not 

being  able  to  supply  the  complete  order  of  echelon  i-1 
in  the  nth  period.  This  covers  the  additional  expected 
cost  resulting  from  echelon  i-1  not  being  able  to  achieve 
its  optimal  level  of  stock  for  the  period.  The  penalty 

A1  is  obtained  from  the  solution  of  echelon  i-1. 
n 

We  now  present  an  example  of  a two- installation  problem.  The 
optimal  policies  are  calculated  for  two  periods. 


3.4 


Two- installation. 
Let  u . 

l 

* , 

hi(ui) 

hi(ui) 

* , 

"AV 

P1<u1) 

Suppose 


Series-configuration  Example 

= inventory  on  hand  in  echelon  i at 
the  end  of  a period. 

= holding  cost  at  installation  i. 

= holding  cost  at  echelon  i. 

= shortage  cost  at  installation  i. 

= shortage  cost  at  echelon  i. 
h1(u1)  = 2.2u1 
h2(u2)  = 2.0u2 


22  - 


T-392 


P*(ux)  = 72ux 

•k 

P2(u2)  = 5u2 
(^(z)  = 5z 

c^(x)  = 50z  + 30 

a = 1 , 

where  all  costs  are  in  dollars. 

Since  this  model  determines  optimal  echelon  policies,  the  period 
cost  , i = 1,2  , must  be  based  upon  echelon  costs,  not  installation 
costs.  The  holding  and  shortage  costs  at  echelon  2 are  merely  these 
respective  costs  at  installation  2.  Thus, 

* 

h2(u2)  = ^2 (u2>  = 2u2  » and 
* 

P2(u2)  = p2(-u2^  = 5u2  ' 

To  find  the  holding  cost  at  echelon  1,  it  is  noted  that  any  units  on  hand 
at  echelon  1 at  the  end  of  a period  have  already  been  charged  $2  holding 
cost  because  they  were  counted  in  the  stock  of  echelon  2 that  was  on  hand 
at  the  end  of  the  same  period.  The  same  reasoning  applies  to  the  penalty 
cost  at  echelon  1.  Thus, 

h^(u^)  = h-^u^)  ~ h^u^)  = 2,2ui  “ 2-Ou^  = 0.2u^  , and 

pi(ui>  = p1(u1)  - p2^ui^  = 72ui  ~ 5ui  = 67ui  • 

If  there  were  N installations  arranged  in  series,  the  holding  and  shortage 
costs  at  echelon  N would  be  obtained  the  same  way  as  for  echelon  2 in 
this  example,  and  the  costs  for  echelons  N-l,  N-2,...,  1 would  be  obtained 

using  the  cost-added  concept  that  was  used  for  echelon  1 in  this  example. 

Note  that  the  cost-added  concept  implies  that  h , > h and  p , > p 

n-l  — n rn-l  — n 

for  n = 2,3, .. .N  . 
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Suppose  the  demand  cj)(t)  for  this  example  follows  a Poisson 
distribution,  and  the  mean  demand  is  one  unit  per  period.  The  period  cost 
is  found  by  evaluating  the  following  expression: 
i 


l,  (y„)  = £ h 

1 n t=0  1 


(y*-t)4>(t)  + I 

t=y1+l 

7n 


p.  * (t-y^WO  » y1  > o , 

l n n — 


i=l,2,  n = 1, 2, . . . 


= £ P,  * (t-y1)(f>(t)  , y1  < 0 , i = 1,2,  n = 1,2,...  . 

t=0  1 


Note  that  the  functions  I.^Cy^)  > i = 1»2  , satisfy  the  assumptions  of 
convexity  because  the  holding  and  shortage  cost  functions  are  convex. 

~1  ~2 

We  now  evaluate  the  recursive  relations  D and  D for  n = 1 

n n 

and  n = 2 , from  which  the  optimal  policies  will  be  obtained.  The  results 
appear  in  Tables  1 and  2,  which  are  read  from  right  to  left.  This  follows 
the  way  dynamic  programming  is  used  to  find  the  optimal  policy  of  the  last 
period  (period  1)  and  then  moves  backwards  in  time,  calculating  successively 
the  optimal  policies  for  periods  2,3,...,n  . Wherever  possible  the  notation 
defined  earlier  is  used  to  head  the  columns  of  these  tables.  Each  table  is 
split  into  two  parts,  the  lower,  representing  echelon  1,  and  the  upper, 
representing  echelon  2. 

Table  1 shows  the  various  possible  costs  for  period  1 at  different 
stock  levels.  The  actual  interpretation  of  stock  level  [Columns  (1)  and 

(8) ],  that  is,  whether  it  is  stock  on  hand  before  or  after  receiving  an 
order,  depends  on  the  cost  column  under  consideration.  Columns  (2)  and 

(9)  of  Table  1 show  zero  cost  because  this  is  the  last  period  of  operation 
and  there  are  no  future  costs.  Successive  differences  of  the  entries  in 
Column  (4)  are  shown  in  Column  (5).  These  are  used  to  locate  the  critical 

,1 


number  S, 


It  was  shown  in  the  description  of  the  classical  single- 


installation model  without  a fixed  cost  of  ordering  that  the  critical 
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Table  1 

PEK10U  1 SOLUTION  OF  TWO-INSTALLATION,  SEKIES-CONFICUKATION  EXAMPLE 


0\ 

Hj(x^)  - 

«l(*) 

Assuming  Always 

. 2 

.2,  2* 

OJ 

-2 , 2 

">1> 

Order  When 

wj(xj-l) 

Total  Expected 
Cost  Before 
Order ing 

W 

A1(xJ) 

a l l>0(y  -t )(|i(t ) 

t-0 

Level 

(15) 

(14) 

(13) 

(12) 

(11) 

(10) 

(9) 

(8) 

12.00 

12.00 

2.00 

2.00 

1.97 

1.87 

1.25 

-12.41 

12.00 

12.00 

0.00 

0.00 

n 

10.00 

10.00 

10.00 

0.00 

0.00 

8.00 

8.00 

8.00 

0.00 

0.00 

■ 

■SSI 

6.03 

6.03 

6.03 

L.OO 

0.00 

H 

mm 

*'♦ . 16 

4.16 

4.16 

0.00 

0.00 

■■ 

itJ 

2.91 

2.91 

2.72 

0.19 

0.00 

MM 

15.32 

55.02 

122.02 

15.32 

55.02 

135.02 

-39.70 
-67 . 00 

15.32 

55.02 

122.02 

2.57 

5.00 

10.00 

12.75 

50.02 

112.02 

0.00 

' 0.00 

0.00 

1 

0 

-1 

185.02 

185.02 

189.02 

15.00 

1 74 . 02 

0.00 

235.02 

235.02 

256.02 

20.00 

236.02 

0.00 

285.02 

285.02 

323.02 

25.00 

298.02 

0.00 

Hi 
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number  S satisfied  W(S  ) = -c  . Since  we  have  shifted  to  discrete 
n n 

valued  functions,  the  method  of  finite  differences  must  be  used  to  find 

S'*"  . It  can  be  shown  that  the  critical  number  S"1  satisfies  the  relation 
n n 

wW)  - Wi(Sn-1)  < ~ C1  < W„<Si+1>  - • <8> 

nn  nn  — i — nn  nn 

Successive  differences  are  shown  in  Column  (5)  until  the  inequality  (8) 
is  satisfied.  In  the  case  of  echelon  1 during  period  1 , 


-5.19  £ -c1  = -5  < -1.08  , 


thus  S^  = 3 . 


Column  (6)  of  Table  1 shows  the  total  expected  cost  after  ordering, 
for  various  stock  levels  before  ordering.  Since  it  is  optimal  to  have  a 

stock  level  of  three  units  at  the  beginning  of  period  one,  if  x^  < 3 , we 

purchase  3 - x^  units  and  add  the  cost  of  this  purchase,  5 • (3-x^)  , to 
~1 

D^(3)  . Column  (7)  shows  the  cost  of  not  being  able  to  bring  the  stock  level 
1 2 

up  to  S^  . This  cost,  A^  , is  equal  to  Column  (4)  minus  Column  (6)  and 

is  the  penalty  charged  to  echelon  2 as  a result  of  installation  2 back- 
ordering installation  1. 

The  echelon  2 section  of  Table  1 shows  the  costs  involved  in  obtaining 
2 2 

the  critical  numbers  (S^,s^).  Much  of  this  section  is  similar  to  the 

echelon  1 section.  Only  the  differences  are  commented  upon.  Column  (10) 

2 

is  equal  to  Column  (7).  In  Column  (12)  the  function  includes  the 

2 2 

additional  cost  A^  . Column  (13)  shows  that  = 0 . Column  (14)  is 

obtained  by  assuming  an  order  is  always  placed  to  bring  the  stock  level  up 
2 

to  zero  if  < 0 . In  comparing  Columns  (14)  and  (12),  note  that  it  is 

2 2 

not  always  wise  to  order  when  x^  < 0 . When  x^  = -1  , the  total  cost 
after  ordering  is  $135.02,  while  if  we  had  not  ordered,  the  total  cost  would 
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have  been  $122.02.  The  stock  level,  s^  , at  which  the  costs  in  Column  (14) 

are  less  than  or  equal  to  the  costs  in  Column  (12),  is  the  level  at  which 

2 

ordering  should  begin.  In  this  case,  s^  = -2  . Column  (15)  shows  the 

total  expected  cost  after  ordering,  for  various  stock  levels  before  order- 
ing. The  entries  of  Column  (15)  are  the  minimum  of  the  entries  in  Columns 
(12)  and  (14). 

Table  2 shows  the  determination  of  the  critical  numbers  for  period  2. 
Columns  (6)  and  (15)  of  Table  1 have  been  brought  back  one  period  in  time 
to  become  Columns  (2)  and  (9),  respectively,  of  Table  2.  The  calculations 
in  Table  2 are  obtained  as  described  for  Table  1. 

The  critical  numbers  for  the  last  two  periods  of  this  problem  are: 
sj  = 3 , (sj.s*)  = (0,-2) 
s2  = 3 . (S^.s*)  = (2’0)  ' 

If  we  are  interested  in  using  this  two-installation,  series-configuration 
example  in  the  design  of  a multi-echelon  system,  which  would  be  in  business 
for  only  two  periods,  and  if  there  are  zero  units  on  hand  at  each  installa- 
tion at  the  beginning  of  period  2 before  an  order  is  placed,  then  a_ 

for  this  echelon  structure  is 

5^(0)  + D^O)  = $23.92  + $165.95  = $139.87  . 

The  calculation  of  the  critical  numbers  for  this  example  were  carried 
out  for  20  periods.  The  critical  numbers  for  the  20th  period  are 

s20  = 5 ’ ^S20’S20^  = (7,;L)  ’ 

The  total  expected  discounted  (inventory)  cost  of  this  echelon  structure 
for  this  product,  that  is,  a^  , when  it  is  assumed  zero  units  are  on  hand 

at  each  installation  before  an  order  is  placed  at  the  beginning  of  the  20th 
period,  is 
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Table  2 

PERIOD  2 SOLUTION  OP  TWO- INSTALLATION , SEKIES-CONP1CUKATION  EXAMPLE 


B2<*2> 

(15) 

“2 

Assuming  Always 
Order  When 

2 < c2 
*2  < 

(14) 

? 7 

W“(x2)  - 

U^xiJ-l) 

(13) 

w^-) 

Total  Expected 
Cost  Before 
Ordering 
(12) 

L2 (y2> 

(11) 

2 2 
A2(*2> 

(10) 

a £ D^(y^-t)iKt) 

t-0 

(9) 

Stock 

Level 

(8) 

2>.02 

22.02 

3.91 
3.53 
2.00 
-2.95 
~ -20.42 
-53.70 

22.02 

12.00 

0.00 

10.02 

7 

L8. 11 

18.11 

18.11 

10.00 

0.00 

8.11 

6 

14.56 

14.58 

14.58 

8.00 

0.00 

6.58 

5 

12.58 

15.53 

35.95 

89.65 

12.58 

15.53 

35.95 

115.95 

12.58 

15.53 

35.95 

89.65 

6.03 

4.16 

2.72 

2.57 

0.00 

0.00 

5.19 

22.74 

6.55 

11.37 

28.04 

64.34 

4 

3 

2 

1 

165.95 

165.95  «i 

189.22 

5.00 

65.02 

119.20 

0 

215.95 

215.95 

322.26 

10.00 

132.02 

180.24 

-1 

265.95 

265.95 

449.04 

15.00 

199.02 

235.02 

-2 

315.95 

315.  95 

571.04 

20.00 

266.02 

285.02 

-3 

365.95 

365.95 

693.04 

25.00 

333.02 

335.02 

-4 

ECHELON  2 


mm 

W*(x*)  - 

w*(-) 

m 

~1  1 

m 

wX-i) 

Total  Expected 
Cost  Before 

a E (yi-t)4>(t ) 

t-0 

Stock 

Level 

mm 

Ordering 

mm 

(6) 

(5) 

(4) 

KB 

(2) 

(i) 

0.00 

0.19 

2.28 

1.20 

1.08 

7 

0.00 

2.11 

1.00 

1.11 

6 

0.00 

HiM 

-0.40 

2.51 

0.84 

1.67 

5 

0.00 

4.28 

-1.77 

4.28 

0.88 

3.40 

4 

0.00 

5.19 

8.92 

13.92 

-4.64 

-10.19 

8.92 

19.11 

1.96 

7.15 

6.96 

11.96 

3 

2 

22.74 

18.92 

41.66 

24.70 

16.96 

1 

65.02 

23.92 

88.94 

66.98 

21.96 

0 

132.02 

28.92 

160.94 

133.98 

26.96 

-1 

199.02 

33.92 

232.94 

200.98 

31.96 

-2 

266.02 

38.92 

304 . 94 

267 . 98 

36.98 

-3 

133.02 

43.92 

376.94 

334 . 98 

41.96 

-4 

ECHELON  1 
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°20(0)  + ^20(0)  = $1438'17  • 

The  series  configuration  of  installations  does  not  represent  a very 
realistic  multi -echel  m system.  If  the  Clark  approach  is  to  be  used  to 
obtain  the  a„  for  the  design  model,  it  will  have  to  be  applied  to  the 

more  realistic  arborescent  configuration  of  installations.  The  next 
section  shows  how  this  approach  is  used  to  find  the  critical  numbers  and 
the  a^j  for  installations  arranged  in  arbjrescence  and  extends  the  method 
to  allow  for  exogenous  demand  at  any  installation. 

3.5  Finding  a„  for  Arborescent-configuration  Structures 

The  purpose  of  this  section  is  to  describe  how  the  Clark  approach 
presented  in  the  previous  sections  is  applied  to  a multi-echelon  situation 
in  which  the  installations  are  arranged  in  an  arborescent  structure. 


Figure  5 shows  the  simplest  multi-echelon,  arborescent-configuration 
system,  where  the  two  installations  at  the  lowest  level  of  the  system,  A^ 

and  A2  , are  fed  by  installation  B . 
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The  situation  illustrated  in  Figure  5 will  be  used  to  describe  how  Clark's 
approach  is  applied  to  arborescent  configurations. 

i 

The  assumptions  in  the  previous  section  for  the  series-configuration 
system  apply  here  with  the  following  exception.  Demand,  exogenous  to  the 
system,  occurs  at  all  the  installations  at  the  lowest  level  of  the  system, 
not  at  only  one  installation.  The  notation  defined  in  the  previous  section 
will  be  used  in  this  section.  A more  cumbersome  method  of  labeling  the 
installations  is  employed  so  that  the  level  as  well  as  the  installation  is 
easy  to  identify. 

Using  Clark's  approach  to  calculate  the  inventory  policies  for  the 
kth  period  for  each  echelon  of  the  system,  we  would  proceed  as  follows. 

(1)  Assume  installations  A^  and  A2  are  independent. 

Determine  the  inventory  policy  for  A^  by  applying 
the  same  method  that  was  used  for  echelon  1 in  the 
series  example,  that  is,  assume  stock  is  available 
at  B to  satisfy  the  order  from  A^  . Also, 

B Ai 

calculate  A^’  1 , the  additional  cost  that  would  be 
experienced  by  installation  A^  if  installation  B 

cannot  satisfy  the  demand  of  A^  during  the  kth  period. 

(2)  Repeat  Step  (1)  for  installation  k^  , that  is, 
calculate  its  inventory  policy  assuming  B can  satisfy 
the  demand  of  A^  . Calculate  A^,A2  . 

(3)  Construct  the  additional  shortage  penalty  A“  from 

B Ai 

the  marginal  per  unit  costs  associated  with  A^’  x 
B Ai  B 

and  A ’ . The  penalty  A is  an  additional  period 

K K 

cost  for  echelon  B during  the  kth  period,  in  the  same 

2 

manner  that  the  penalty  A^  applied  to  echelon  2 of 
the  example  in  the  previous  section. 
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(4)  Determine  the  inventory  policy  for  echelon  B.  Since 

a fixed  cost  of  ordering  is  allowed  at  installation  B, 
the  highest  installation  in  this  system,  the  method 

B B 

for  determining  the  two  critical  numbers  (S^,s^) 

follows  the  method  used  for  echelon  2 of  the  series- 
configuration  example. 


In  short,  the  method  outlined  above  says  to  calculate  the  inventory 
policies  for  each  echelon  separately  by  the  methods  used  for  the  classical 
single- installation  situation.  This  procedure  can  be  extended  easily  to 
apply  to  a situation  where  B feeds  more  than  two  installations  or  to  a 
system  containing  more  than  two  levels. 

Recall  that  C (x  1,  x ^ , x ) is  the  minimum  total  expected 
n n n n 

discounted  cost  for  the  entire  system  of  this  example  at  the  beginning  of 
period  n before  an  order  is  placed.  If  the  inventory  policies  calculated 
by  this  method  are  to  minimize  the  total  expected  discounted  cost  of  the 
system,  it  is  necessary  that 


~ , Ai  At  B.  ~Ai  Ai  . . £A9.  A?x  , ~B.  B. 

C (x  x z,  x ) = D X(x  ■*■)  + D ^(x  ^)  + D (x  ) . 


Unfortunately,  it  has  been  shown  by  Clark  and  Scarf  (1960)  that  decomposition 
(9)  does  not  always  hold  for  the  arborescent  configuration  of  installations. 
They  give  an  argument  for  such  a decomposition  by  assuming  that  the  stock 
levels  at  the  lowest  installations,  in  this  case  A^  and  k^  , are  not 

out  of  balance.  By  this  they  mean  that  for  a given  period  the  ratios  of 
stock  on  hand  after  ordering  to  expected  demand  are  approximately  the  same 
for  A^  and  A^  . In  practical  situations  they  feel  that  the  stock  levels 

at  the  lowest  installations  are  rarely  out  of  balance,  and  therefore  they 
conclude  that  Clark's  approach  gives  excellent  approximations,  if  not 
optimal  solutions,  for  the  best  inventory  policies  of  an  arborescent 
configuration. ^ 

* A /\ 

The  notations  C and  D will  continue  to  be  used  in  the  remainder  of 
this  section,  even  though  it  cannot  be  proven  that  these  are  the  minimum  costs. 


L 
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If  the  holding,  shortage,  and  transportation  costs  per  unit  at 
each  of  the  lowest  level  installations  are  of  the  same  order  of  magnitude, 
then  it  is  felt  that  the  stock  levels  at  the  lowest  installations  generally 
will  be  in  balance.  Experience  with  Clark’s  approach  for  obtaining  inven- 
tory policies  of  installations  arranged  in  an  arborescent  configuration 
indicates  that  his  method  is  probably  quite  good.^ 

An  assumption  was  made  which  should  not  be  overlooked.  It  was 
assumed  in  Step  (1)  that  installations  and  were  independent. 

This  means  that  transshipments  between  A^  and  are  not  allowed.  Thus, 

an  installation  will  be  allowed  to  receive  stock  only  from  its  designated 
supplier  at  a level  higher  in  the  system. 

Clark’s  approach  is  now  used  to  obtain  the  inventory  policies  for 
a three- installation,  arborescent-configuration  problem. 


3.6  Three-installation,  Arborescent-configuration  Example 
Suppose 


(UA  > 
A1 


(UA  > 
2 


VV 


(UA  > 
A1 


(UA  } 
2 


PB(UB) 


2.2u. 


2.1u. 


2 . Ou 


B 


72u. 


69u 


5 . Ou 


B 


1 


2 


^Clark  (1960)  used  this  approach  to  obtain  the  inventory  policies 
that  were  used  in  the  simulation  of  a large  arborescent-configuration 
system.  In  personal  conversations  with  him  he  indicated  that  the  policies 
produced  for  this  simulation  were  excellent. 


- 32  - 


T-392 


c (z)  = 5z 
A1 


c (z)  = 3z 
A2 


c (z)  = 50z  + 30 

D 


a = 1 


Using  the  cost-added  concept  described  in  the  series-configuration  example, 
the  echelon  holding  and  shortage  costs  are: 

hA  (uA  } = °‘2uA 
1 1 1 


hA  (uA  } = °,1UA 
2 2 2 


W 


= 2.  Ou 


B 


PA  (UA  > 
A1  A1 


= 67u , 


PA  (uA  } 
A2  2 


= 64u . 


w 


5uB  * 


Suppose  the  demands  at  installations  A^  and  A^  each  follow  a Poisson 

distribution,  with  mean  demand  one  unit  per  period.  The  demand  experienced 
by  echelon  B is  the  convolution  of  the  demands  experienced  by  the  two 
installations  fed  by  B.  Therefore,  the  demand  at  echelon  B follows  a Poisson 
distribution,  with  a mean  demand  of  two  units  per  period. 

The  calculations  for  the  inventory  policies  at  each  echelon  for 
period  1 are  now  described.  We  follow  the  four-step  procedure  given  above. 
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(1)  The  costs  and  demand  distribution  for  echelo;.  are 

identical  to  those  given  for  echelon  1 of  the  series- 
configuration  example.  Therefore,  the  echelon  1 section 
of  Table  1,  Columns  (1)  - (7),  gives  the  calculations  for 

period  1 of  echelon  in  this  example  and  SA^  = 3 . 

Column  (7)  of  Table  1 is  . 

(2)  Table  3 gives  the  calculations  for  period  1 of  echelon 
A^  . The  arrow  indicates  the  location  of  the  critical 

number  and  SA^  = 3 . 

(3)  We  now  show  how  A^’A1  and  A^’A^  are  used  to  construct 

g 

A^  . If  installation  B cannot  satisfy  the  demand  from  A^ 
and  A^  , it  has  a choice  of  backordering  A^  , A^  , or 

both  A^  and  A^  • It  desires  a plan  for  backordering 

that  will  penalize  it  the  least.  For  example,  if  B is 
short  one  unit  in  the  first  period,  it  is  clear  by  examining 
Column  (7)  in  Table  1 and  Column  (7)  in  Table  3 that  it  is 
cheaper  to  backorder  A^  one  unit,  where  the  penalty  is 
$0.19,  than  to  backorder  A^  one  unit,  where  the  penalty 

is  $2.04.  Similarly,  if  B were  short  three  units  in 
period  1,  it  has  the  following  four  alternatives: 

(a)  backorder  A^  three  units  at  a cost  of 
$50.02 

(b)  backorder  A^  two  units  and  A^  one  unit 
at  a cost  of  $12.75  + $2.04  = $14.79 

(c)  backorder  A^  one  unit  and  A^  two  units 
at  a cost  of  $0.19  + $15.87  = $16.06,  or 
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(d)  backorder  three  units  at  a cost  of 


$53.29. 

Thus,  B would  choose  alternative  (b)  in  order  to  minimize 
the  cost  of  not  being  able  to  bring  the  lower  installations 
up  to  their  desired  stock  levels.  It  is  assumed  that  instal- 
lation B will  use  the  minimum  cost  backordering  alternative, 
and  furthermore,  that  B generally  has  a choice  regarding  how 
much  it  will  backorder  each  lower-level  installation.  To 

g 

determine  the  penalty  A^  for  all  possible  backordering  situa- 
tions, we  first  calculate  the  marginal  per  unit  costs  associa- 
ted with  A^,A1  and  A^’A2  } called  marginal  A^,A1  and  mar- 
ginal A^,A2  in  Table  4.  These  give  the  additional  cost  for 
backordering  one  more  unit  at  and  , respectively. 

The  marginal  penalties  are  then  ranked  in  value,  starting  with 
the  smallest.  The  ranked  penalties  are  then  successively 

g 

added  to  form  A^  . This  procedure  guarantees  that  the  least 

cost  combination  in  backordering  the  lower-level  installation 
is  used.  These  calculations  are  shown  in  Table  4.  The  posi- 

BA  BA 

five  ranked  marginal  A^’  1 and  A^’  2 start  at  stock  level  5 
because  if  echelon  B has  less  than  six  units,  then  either 
or  A2  are  below  their  desired  level  of  three  units  each. 

g 

(4)  Using  the  A“(*)  from  Table  4,  the  inventory  policy  for 

echelon  B at  the  beginning  of  period  1 can  be  determined. 

The  calculations  are  shown  in  Table  5 and  follow  the  calcula- 
tions for  echelon  2 of  the  series-configuration  example,  given 
in  Table  1.  The  critical  numbers  as  indicated  by  the  arrows 

are  (S®,s®)  = (0,-2)  . 
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If  interest  was  only  in  one  period  of  operation  for  this  three- 
installation  arborescent-configuration  example,  and  if  the  on-hand  inventory 
at  each  installation  at  the  beginning  of  period  1 before  ordering  is  zero 
units,  then  the  total  expected  discounted  cost  of  the  system  would  be 


D^1(0)  + D^2(0)  + dJ(0)  = $16.96  + $10.69  + $113.32  = $140.97 


This  would  be  the  inventory  cost  a_  for  the  design  model.  The  calcula- 
tions given  here  for  one  period  were  carried  out  for  20  periods.  The  cri- 
tical numbers  for  the  20th  period,  that  is,  with  20  periods  of  business 
remaining,  are 

S20  = 5 * S20  = 5 * (S20’S20}  = (11’3)  * 


The  total  expected  discounted  cost  for  20  periods  is  $2681.29,  and  this 
would  represent  the  inventory  cost,  a^  , in  the  design  model  if  we  were 

designing  a system  to  last  for  20  time  periods. 


3.7  Shadow  Installations 


■ 


The  method  just  described  for  obtaining  near  optimal  inventory 
policies  for  N installations  having  an  arborescent  configuration  will  be 
satisfactory  for  many  of  the  echelon  structures  to  be  considered  by  the 
design  model.  But  there  is  one  situation  which  cannot  be  handled  directly 
by  the  Clark  approach  for  arborescent  configurations.  This  is  the  situation 
in  which  one  or  more  of  the  lowest-level  installations,  where  exogenous 
demand  exists,  is  not  included  in  the  echelon  structure.  This  would  be  the 
situation  if  one  or  both  of  installations  or  were  removed  from 

the  system  of  Figure  5.  For  example,  Figure  6 shows  this  system  with 
installation  A^  removed.  Although  installation  A ^ is  not  physically 

included  in  the  echelon  structure  of  Figure  6,  it  is  represented  on  the 
diagram  by  dashed  lines  because  demand  still  exists  where  this  installa- 
tion would  be  located  if  it  were  included  in  the  structure.  We  call 
installation  A2  a shadow  installation.  Any  time  an  echelon  structure 
does  not  include  one  of  the  lowest-level  installations,  that  installation 
is  called  a shadow  installation. 
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In  the  echelon  structure  of  Figure  6,  demand  must  be  satisfied 

by  installation  B on  an  individual  basis,  that  is,  it  is  a "mail-order" 
demand. ® 

Source 

1 


t i 

»X  D2 

Figure  6.  Arborescent  configuration  with 
a shadow  installation. 

Suppose  we  desire  to  determine  the  inventory  policies  for  echelons 

and  B , taking  into  account  that  mail-order  demand  exists  at  the 

location  where  installation  A^  would  be  placed  if  it  were  included  in  the 

system.  To  do  this,  we  need  a method  which  includes  the  transportation 
cost  from  installation  B to  the  customers  represented  by  , and  a penalty 

cost  at  B for  not  satisfying  the  demand  at  . There  is  no  need  for  a 

holding  cost  at  since  this  installation  is  not  in  the  system.  Let 

— -k 

p.  = the  shortage  cost  per  unit  for  customers  at 
the  shadow  installation  i. 


c^  = the  transportation  cost  per  unit  in  getting 
an  item  to  the  customers  of  shadow  installa- 
tion i from  their  supplier. 


We  use  the  phrase  "mail-order"  in  referring  to  these  individual 
demands,  but  the  manner  in  which  they  are  processed  is  completely  arbitrary, 
as  long  as  they  are  processed  on  an  individual  basis.  Therefore,  if  the 


T-392 


The  echelon  structure  of  Figure  6 is  used  to  describe  the  method  for 
obtaining  the  inventory  policies  and  a^  for  a multi-echelon  configuration 

with  a shadow  installation.  The  procedure  for  the  kth  period  is: 

(1)  Assume  demand  can  only  be  satisfied  by 

installation  B,  that  is,  it  cannot  be  satisfied 
by  installation  . Determine  the  inventory 

policy  for  A^  and  A®’  1 using  Clark’s  approach 

for  arborescent  configurations. 


(2) 


(3) 


No  calculations  are  required  for  echelon  A ^ , 
since  installation  A^  does  not  exist. 


To  construct  A^  , it  is  necessary  to  consider  a 

penalty  at  echelon  B for  not  being  able  to  satisfy 
the  demand  , as  well  as  the  demand  from  A^  . 

Thus,  we  need  something  similar  to  the  marginal  per 


unit  cost  of  A! 


b,a2 


Whenever  a customer  from 


demand  is  backordered  by  installation  B,  he 

* 

incurs  a cost  of  pA  . But  part  of  this  cost  has 

a2 

* 

already  been  applied  by  p ; therefore,  using  the 
cost-added  concept,  we  let  p^  equal  the  cost  per 

k 

unit  added  by  penalty  p . Then,  marginal 

A2 

. B , Ao  — 

A * = D 


(4)  The  inventory  policy  for  echelon  B is  determined  in 

the  same  way  that  it  was  found  in  the  arborescent 

configuration,  with  the  following  exception.  The 

period  cost,  !,„(•)  , must  include  the  additional  cost 
B 

of  meeting  the  expected  demand  of  D^  , , or  in 

this  case 
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y£ 

= "A2A2  + tlQ 

+ ^ PB(t”yk)<Vt)  » ^k  - ° ’ k = 

t=yk+l 

= CA  ^2  + Z pb^t-yk^B^^  , < ° , k = 1,2 

This  procedure  is  now  used  to  obtain  the  inventory  policies  for  a 
numerical  example  of  the  arborescent  configuration  with  a shadow  installation 
shown  in  Figure  6. 


3.8  Shadow  Installation  Example 
Suppose 

y V - 2-y 

hB<v  ■ 2-% 

Yy  ■ 72UA! 
v - 78 


W ‘ 5uB 


c (z)  = 5z 
A1 


c.  = 10 
A2 


c_(z)  = 50z  + 30 

b 


a = 1 . 
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Using  the  cost-added  concept,  the  echelon  holding  and  shortage  costs  are: 

hA  (UA  } = °*2uA 
1 1 1 

W = 2 ’ °UB 

= 67uax 

PA  " 73 
2 

PB(UB)  = 5uB  • 

Suppose  the  demands  and  are  Poisson,  with  mean  demand  one  unit 

per  period. 

The  calculations  of  the  inventory  policies  at  echelons  A^  and  B 

for  period  1 are  now  described,  following  the  procedure  given  for  multi- 
echelon configurations  with  a shadow  installation. 

(1)  The  costs  and  demand  distribution  for  echelon  A^ 

are  identical  to  those  given  for  echelon  1 of  the 
series-configuration  example.  Therefore,  the  eche- 
lon 1 section  of  Table  1,  Columns  (1)  - (7),  gives 
the  calculations  for  period  1 of  echelon  A^  in 
Ai 

this  example  and  = 3 . Column  (7)  of  Table  1 

is  A®’Al  . 

(2)  Installation  A ^ is  not  included  in  the  configuration, 

A? 

therefore  Sn ^ is  not  calculated. 

B 

(3)  A^  is  constructed  for  this  example  using  the  same 

method  applied  in  the  three- installation  arborescent- 
configuration  example,  with  the  exception  that  the 
B Ao  - 

marginal  A..  ’ *■  = p = 73  , for  all  stock  levels.  The 
i a2 
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B i 

ranking  of  the  marginal  A ’ , i = A^  and  , and  the 

JJ 

construction  of  A^  are  shown  in  Table  6.  The  positive 
B A9 

marginal  A ’ 2 begin  at  stock  level  0 because  we  expect 

a mail-order  demand  from  D ^ of  one  unit  per  period  and, 

therefore,  want  to  penalize  echelon  B for  holding  less  than 

B Ai 

one  unit  of  stock.  The  positive  ranked  marginal  A^ ’ 1 
B A 9 

and  A^’  z start  at  stock  level  3 because  if  echelon  B 

has  less  than  four  units,  then  either  A^  will  be  below  its 

desired  stock  level  of  three  units  or  the  expected  demand  of 
one  unit  at  D£  will  not  be  met. 

g 

(4)  Using  the  A^(*)  from  Table  6,  the  inventory  policy  for 

echelon  B at  the  beginning  of  period  1 can  be  determined. 

The  one  exception  to  the  method  used  for  the  arborescent- 
configuration  example  is  that  L (•)  must  include  the  addi- 

D 

tional  cost  of  meeting  the  expected  demand  of  for  each 

value  of  stock  level.  This  cost  equals  the  product  of  the 

expected  demand  and  the  transportation  cost,  c , or  in  this 

2 

case  $10.  The  costs  and  distribution  for  demand  at  echelon  B 
in  this  example  are  identical  to  those  used  for  echelon  B 
of  the  three-installation  arborescent-configuration  example. 
Therefore  LD(*)  for  this  example  is  the  L (•)  found  in 

D D 

Table  5 plus  $10.  The  calculations  for  echelon  B follow  the 
calculations  for  echelon  2 of  the  series  configuration  example, 

B B 

given  in  Table  1.  The  critical  numbers  are  (S^,s^)  = (1,-1). 

The  determination  of  the  inventory  policies  for  periods  2,3,...  , 
follow  the  method  described  for  period  1.  The  calculations  for  this  example 
have  been  carried  out  20  periods,  and  the  critical  numbers  for  the  20th 
period  are: 
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Table  6 


,B 


DETERMINATION  OF  A“  FOR  ARBORESCENT-CONFIGURATION 


EXAMPLE  WITH  A SHADOW  INSTALLATION 


Stock 

Level 

ab’ai 

1 

Marginal 

ab,ax 

Marginal 

ab,A2 

Ranked  Marginal 

B,Ai  , B ,A? 

A^  1 and  A^  z 

“l 

7 

0.00 

0.00 

0.00 

0.00 

0.00 

6 

0.00 

0.00 

0.00 

0.00 

0.00 

5 

0.00 

0.00 

0.00 

0.00 

0.00 

4 

0.00 

0.00 

0.00 

0.00 

0.00 

3 

0.00 

0.00 

0.00 

0.19 

0.19 

2 

. 0.19 

0.19 

0.00 

12.56 

12.75 

1 

12.75 

12.56 

0.00 

37.28 

50.02 

0 

50.02 

37.28 

73.00 

62.00 

112.02 

-1 

112.02 

62.00 

73.00 

62.00 

174.02 

-2 

174.02 

62.00 

73.00 

62.00 

236.02 

-3 

236.02 

62.00 

73.00 

62.00 

298.02 

-4 

298.02 

62.00 

73.00 

62.00 

360.02 

r 
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1 


S20  ' 5 “d  <S5o-‘5o>  ' <9-2)  ' 

If  there  are  zero  units  on  hand  at  each  installation  before  ordering  in 
the  20th  period,  the  total  expected  discounted  cost  for  20  periods  is 
$2708.11.  This  would  represent  the  inventory  cost,  a^  , in  the  design 

model  if  we  were  designing  a distribution  system  to  last  for  20  periods 
and  were  interested  in  this  echelon  structure.  Note  that  the  total  expected 
discounted  cost  of  operating  the  three-installation  arborescent-configuration 
was  $2681.29.  Therefore,  in  comparing  these  two  echelon  structures,  we  might 
conclude  that  the  mail-order  operation  is  a poor  way  to  do  business  for  this 
situation.  3ut  this  does  not  take  into  account  the  inventory  costs  for  other 
products  or  the  fixed  cost  of  an  installation. 


4.  THE  DESIGN  PROBLEM  WITH 
STORAGE  SPACE  CONSTRAINTS 


In  this  section  we  generalize  the  design  problem  formulated  in 
Section  2 to  allow  for  storage  space  constraints  at  the  various  facilities. 
The  resulting  model  yields  a 0-1  linear  programming  problem,  which  includes 
the  previously  discussed  design  problem  without  capacity  constraints  as  a 
special  case. 

Formulation  of  the  design  problem  with  capacity  constraints  requires 
notation  in  addition  to  that  given  in  Section  2,  ana,  for  convenience,  we 
give  here  all  notation  used  in  the  model.  The  echelon  structures  are  indexed 
by  i,i=l,...,m;  the  products  are  indexed  by  j,j=l,...,n;  the  installations 
are  indexed  by  k,k=l,...,p.  Let 


a 


ij 


the  inventory  cost  of  product  j 
using  echelon  structure  i, 


b^  = the  facility  cost  of  installation  k , 


r^  = the  storage  space  available  at  installation  k , 


d 


ijk 


the  storage  space  required  at  installation  k for 
product  j when  product  j uses  echelon  structure  i. 
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The  d, values  could  be  measured  on  the  basis  of  either  the  maximum 
ljk  


space  required  (over  time)  or  the  average  space  required  (over  time).  The 
former  seems  more  appropriate  for  this  type  of  design  problem,  and  the 
r^  values  must  then  be  determined  accordingly.  Note  that  the 


d..,,  k =l,...,p  , values  are  obtained  as  outputs  of  the  multi-echelon 
lj  K. 


inventory  problem  solved  in  calculating  a_  . Thus  d„^  = ® 


j=l,...,n  if  echelon  structure  i does  not  use  installation  k. 


The  decision  variables  are: 


x„  =1  if  product  j uses  echelon  structure  i , 


0 otherwise. 


y.  = 1 if  installation  k is  used, 
k 


= 0 otherwise. 

The  problem,  henceforth  called  problem  (P) , is  then  to  find  x^.  and  y^ 


values  that 


m n 


minimize  Z Z a. .x. . + Z b,  y, 
i=!  j=l  « ^ k=l  k k 


(10) 


subject  to 


m 


i=l 


m n 


« 


(ID 

(12) 


x..,y  = 0 or  1 for  all  i,j,k  . (13) 

ij  k 


Comparison  of  this  formulation  with  that  of  Section  2 shows  that  the  objective 
functions  (1)  and  (10)  are  equivalent  and  constraints  (2)  and  (11)  are  the 
same.  The  constraints  (12)  serve  to  impose  the  storage  space  limitations  at 
the  same  time  that  they  force  the  use  of  an  installation  if  any  echelon 

must 


structure  using  it  is  chosen.  If  y^  = 0 and  d„^  > ® » the*1 


be  0 in  order  to  satisfy  (12).  Note  that  the  previous  formulation  without 
storage  space  constraints  can  be  easily  obtained  as  a special  case  of 
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I 


(10)  - (13)  by  letting  d..  equal  1 if  echelon  structure  i uses  installation 

1 J K 

k (and  0 otherwise)  and  setting  all  r^  equal  to  n. 

The  design  problem  with  storage  space  constraints  has  been  formulated 
in  this  section  as  a 0-1  LP  problem.  In  the  next  section  we  illustrate 
its  use  for  designing  a multi-product,  multi-echelon  inventory  distribution 
system. 

5.  ILLUSTRATION  OF  MODE’S 

This  section  presents  illustrative  computational  results  using  the 
three  models  given  in  this  paper — the  dynamic  programming  model  to  find 
optimal  inventory  policies  for  a multi-echelon  distribution  system,  the 
design  model  without  storage  space  constraints,  and  the  design  model  with 
storage  space  constraints.  The  section  concludes  with  a discussion  of  some 
problems  yet  to  be  resolved. 

5.1  Illustration 

In  this  illustration  five  alternative  echelon  structures,  which  include 
up  to  three  echelons  and  eight  facilities,  are  being  corsidered  for  the 
design  of  a four-product,  multi-echelon  inventory  distribution  system. 

The  five  echelon  structures  are  shown  in  Figure  7.  Echelon  structure  1 
is  the  full  system,  that  is,  it  consists  of  the  maximum  number  of  facilities 
under  consideration  (a  central  warehouse,  two  regional  warehouses,  and 
five  retail  stores)  arranged  in  three  echelons.  Although  eight  facilities 

Q 

generate  2 = 256  possible  echelon  structures,  only  the  four  additional 

echelon  structures  shown  are  being  considered  for  this  system.  Echelon 
structure  2 is  a one-echelon  system  with  mail  order  delivery  from  the 
central  warehouse  to  customers  located  where  the  retail  stores  would 
be  if  they  were  in  the  system.  Echelon  structure  3 is  a two-echelon 
system  without  a central  warehouse.  In  this  case  the  two  regional  warehouses 
are  supplied  directly  by  the  sources  of  production.  Echelon  structure  4 
is  a single  echelon  system  where  the  retail  stores  are  supplied  with 
products  directly  from  the  sources  of  production,  while  echelon  structure  5 
is  similar  to  echelon  structure  1,  with  the  exception  that  customers  at 
retail  stores  1,  2,  and  3 are  supplied  by  mail  order  from  regional 
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warehouse  6. 

The  dynamic  programming  approach,  described  in  Section  3,  was 
used  to  find  the  multi-echelon  inventory  policies  and  the  resulting 
a„'s  for  the  four  products  of  this  example.  Table  7 summarizes  the 

input  data  needed  to  obtain  these  policies.  It  is  assumed  that  the 
stock  levels  of  this  inventory  system  are  reviewed  quarterly  and  the 
demand  for  each  product  has  a Poisson  distribution,  with  the  mean 
demand  per  quarter  assumed  to  be  the  same  at  each  retail  store.  If 
echelon  structure  1 is  used  to  stock  a product,  it  is  generally  assumed 
that  the  cost  to  transport  the  product  from  the  central  warehouse 
(facility  8)  to  one  of  the  retail  stores  is  6%  of  the  procurement  cost, 
that  is,  the  cost  to  get  the  product  to  facility  8.^  Referred  to  in 
Table  7 as  "normal"  transportation  cost,  this  cost  is  divided  evenly 
between  the  cost  to  move  the  product  between  facility  8 and  one  of  the 
regional  warehouses  (facility  6 or  7)  and  the  cost  to  move  the  product 
between  the  regional  warehouse  and  one  of  the  retail  stores. 

In  the  case  of  echelon  structures  with  shadow  installations,  the 
transportation  cost  to  satisfy  mail-order  demands  will  generally  increase 
significantly.  Table  7 shows  these  costs  as  a percentage  of  the  normal 
transportation  costs  described  above.  Also  shown  are  the  changes  in 
transportation  cost  when  either  the  central  warehouse  or  regional  warehouse 
are  not  used  to  stock  a given  product.  For  products  1 and  2 we  assume 
under  such  a two-echelon  system  that  there  would  be  no  change  from  the 
normal  transportation  cost,  whereas  with  products  3 and  4 it  is  assumed 
that  the  transportation  costs  double  (200%  of  normal  cost)  under  these 
circumstances.  Table  7 also  gives  the  holding  and  stockout  penalty  costs 
used  in  this  illustration.  The  holding  costs  per  quarter  range  from  2% 
to  8%  of  the  total  cost  to  get  the  product  to  the  facility  where  this 
holding  cost  is  incurred.  The  shortage  costs  at  the  retail  stores 


9 

In  two  cases,  for  products  2 and  4 transported  to  retail  stores  4 and  5, 
this  cost  is  assumed  to  be  3%  of  the  procurement  cost. 
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(or  shadow  installations)  vary  from  110%  to  190%  of  the  cost  to  get  a 
product  to  these  stores.  At  facilities  other  than  the  retail  stores,  a 
shortage  cost  equal  to  10%  of  the  facility's  holding  cost  is  applied  to 
cover  the  nuisance  created  by  backordering. 

Finally,  the  following  data  or  characteristics  were  assumed  to  be 
common  to  all  the  products: 


1.  The  life  of  the  system  is  20  quarters  (five  years), 
that  is,  each  product  is  stocked  for  20  quarters, 

2.  there  is  no  inventory  on  hand  at  the  beginning  of 
the  first  period  of  operation, 

3.  the  interest  rate  is  1-1/2%  per  quarter,  resulting 
in  a discount  factor  of  0.985,  and 


4.  the  fixed  cost  of  ordering  is  $100. 

Table  8 gives  the  dynamic  programming  results.  The  a„  represent 

the  total  expected  discounted  cost  for  stocking  product  j under  echelon 
structure  i for  20  quarters,  and  the  inventory  policies  shown  are  for 
the  20th  quarter.  The  computer  program  written  to  generate  these  results 
also  gives  the  order  policies  at  each  facility  for  each  of  the  20  periods. 

It  is  interesting  to  note  that  steady  state  policies  were  generally  reached 
quickly.  With  the  exception  of  three  of  the  20  results,  the  policy  shown 
for  the  20th  quarter  was  reached  by  the  fifth  quarter,  and  sometimes  sooner. 


Two  types  of  inventory  policies,  S and  (S,s),  are  shown  in  Table  8. 
At  the  highest  level  facility  in  a system  (the  facility  that  orders  the 
product  from  the  source  of  production)  it  is  assumed  that  a fixed  cost  of 
ordering  is  incurred,  hence  (S,s)  type  policies  result.  Since  there  is 
no  fixed  cost  of  ordering  at  other  facilities,  S type  policies  are  found. 
The  time  to  calculate  these  critical  numbers  and  to  find  the  inventory 
costs,  a„  , averaged  30  seconds  per  product  for  a given  echelon  structure 
on  an  IBM  360/50  computer. 


Next  we  determine  the  best  design  for  this  four-product  system, 
using  the  model  given  in  Section  2,  that  is,  we  assume  there  are  no 
constraints  on  the  amount  of  space  available  at  each  of  the  facilities 
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where  products  might  be  stored.  The  inventory  costs,  a„  , given  in 

Table  8,  and  the  facility  costs,  b ^ > for  this  problem  are  arrayed  by 

echelon  structure  i in  Figure  8.  The  facility  costs  are  consistent 
with  values  found  in  previous  research  [see  Pinkus  (1971)]  and  take 
into  account  the  fact  that  the  fixed  cost  of  operating  a warehouse,  in 
proportion  to  the  value  of  goods  that  might  be  stocked  there,  is  less 
than  the  fixed  cost  of  operating  a retail  store,  due  to  economies  of 
scale  and  lower  property  values. 


The  solution  to  this  problem  is  to  stock  all  four  products  under 
echelon  structure  5,  which  means  that  each  product  will  be  stocked  at 
retail  stores  4 and  5 and  at  warehouses  6,  7,  and  8.  The  total  (5-year) 
expected  discounted  cost  of  this  solution  is  $766,000.  It  is  interesting 
to  note  that  of  the  four  products,  only  products  1 and  3 use  the  echelon 
structure  that  minimizes  their  inventory  cost.  If  products  2 and  4 were 
to  use  echelon  structure  1,  which  minimizes  their  inventory  cost,  the 
fixed  cost  of  the  system  would  increase  because  products  2 and  4 would 
now  be  stocked  at  retail  stores  1 - 3,  as  well  as  at  facilities  4-8. 
Under  this  solution  the  savings  in  inventory  cost  do  not  compensate  for 
the  increased  facility  cost,  and  the  total  cost  of  the  system  is  $790,000. 

The  solution  to  this  problem  was  obtained  in  11  seconds  on  an 
IBM  370/148  using  the  branch-and-bound  algorithm  described  in  Pinkus, 
Gross,  and  Soland  (1973).  Although  initially  written  to  handle  10 
echelon  structures,  30  products,  and  30  facilities,  the  computer  program 
for  this  algorithm  can  easily  be  expanded  to  handle  larger  values  of  m , 
n , and  p . 


Finally,  we  have  solved  the  same  problem  with  storage  space 
constraints  using  the  model  presented  in  Section  4.  The  storage  space 
available  at  each  facility,  r^  , is  given  in  Table  9,  and  the  amount 


of  storage  space  required  for  each  product  is  given  in  Table  10. 


To 


find 


d , we  multiply  the  per  unit  storage  space  required  to  stock 


product  j (see  Table  10)  by  the  maximum  number  of  units  (S)  that  will 
be  stocked  at  facility  k when  echelon  structure  i is  used  (see  Table  8). 
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Table  9 


STORAGE  SPACE  AVAILABLE 
AT  EACH  FACILITY 


Facility 


(sq.  feet) 


Table  10 

STORAGE  SPACE  REQUIRED  BY 
EACH  PRODUCT 


Product 


Storage  Space 
Required  Per  Unit 
(sq.  feet) 
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For  example,  d.  n , is  800  (that  is,  400  x 2),  d_  „ , is  900  (that  is, 

1,1)1  l)ZtO 

300  x 3),  and  d.  < 0 is  3000  (that  is,  300  x 10). 

D ) 4 , o 

The  solution  to  the  problem  with  space  constraints  is  to  stock 
products  1,  2,  and  4 under  echelon  structure  4,  and  product  3 under 
echelon  structure  2,  that  is,  products  1,  2,  and  4 are  stocked  only 
at  the  retail  stores,  which  are  supplied  directly  by  the  sources  of 
production,  and  product  3 is  stocked  only  at  the  central  warehouse 
(facility  8),  from  which  customers  receive  product  3 by  mail  order. 

The  total  (5-year)  expected  discounted  cost  for  this  solution  is  $779,000. 

In  this  case,  none  of  the  products  are  stocked  under  the  echelon  structure 
that  minimizes  the  inventory  cost  for  a given  product.  This  solution  was 
obtained  in  145  seconds  using  a 0-1  linear  programming  computer  code  run 
on  an  HP  3000. 

5.2  Discussion 

We  conclude  this  paper  with  a discussion  of  some  problems  yet  to 
be  resolved. 

The  design  model  with  storage  space  constraints,  problem  (P) , has 
mn  + p 0-1  variables  and  n + p constraints,  so  the  problem  dimensions 
may  be  fairly  high  for  problems  of  practical  size.  For  example,  with 
m = n = 30  and  p = 20  , problem  (P)  has  920  variables  and  50  constraints. 
The  0-1  LP  computer  code  used  to  solve  the  design  problem  with  storage 
space  constraints  in  the  previous  section  can  handle  only  40  variables 
and  20  constraints.  This  fact,  together  with  the  special  structure  of 
problem  (P) , suggests  that  a specialized  algorithm  for  its  solution  would 
be  much  more  efficient  for  practical  problems  than  the  general  integer 
linear  programming  algorithm  that  was  used  to  get  the  results  in  Section  5.1. 
Such  a special  algorithm  is  presently  being  developed  and  tested. 

There  is  another  aspect  of  this  work  that  needs  to  be  examined.  It 
can  be  shown  that  the  model  with  space  constraints,  problem  (P) , might  not 
obtain  the  best  solution  to  this  design  problem.  This  is  because  we  have 
restricted  ourselves  to  using  the  optimal  inventory  policy  and  resulting 
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a.,  for  a given  product  stocked  under  a given  echelon  structure.  With 


the  introduction  of  space  constraints,  a situation  can  occur  where  by 
making  a reduction  in  the  optimal  value  of  S for  a given  product,  we 
might  be  able  to  satisfy  a space  constraint  that  would  not  otherwise  be 
satisfied.  Such  a reduction  in  S would  yield  a sub-optimal  (greater) 
value  for  a„  , but  this  increased  inventory  cost  could  be  compensated 

for  by  the  reduction  in  facility  cost  resulting  from  enabling  the  product 
to  be  squeezed  into  a set  of  facilities  used  by  other  products. 

There  is  also  the  possibility  of  letting  the  decision  variables, 

x. . , take  on  fractional  values.  A fractional  x..  means  that  part  of 
ij  ij 

product  j is  stocked  under  one  echelon  structure,  and  the  remainder 

under  one  or  more  other  echelon  structures.  It  is  easy  to  see  that  the 

x„  would  always  be  0 or  1 for  the  design  problem  without  storage  space 

constraints,  but  it  is  not  clear  that  this  should  be  the  case  for  problem  (P). 
We  are  in  the  process  of  allowing  for  the  use  of  fractional  x„  , for  the 

space  constrained  design  problem,  and  plan  to  investigate  several  heuristic 
solution  methods  in  order  to  allow  for  the  possiblity  of  using  sub-optimal 
inventory  policies  to  reduce  the  total  cost  of  the  system. 
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